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application to a phase qubit 



N. Didier^ and F. W. J. Hekking^ 

^ Scuola Normale Superiore, NEST and Istituto Nanoscienze-CNR, Pisa, Italy 
^Laboratoire de Physique et de Modelisation des Milieux Condenses, 
Universite Joseph Fourier and CNRS, BP 166, 38042 Grenoble, France 

Macroscopic quantum tunneling of the phase is a fundamental phenomenon in the quantum dy- 
namics of superconducting nanocircuits. The tunneling rate can be controlled in such circuits, where 
the potential landscape for the phase can be tuned with different external bias parameters. Precise 
theoretical knowledge of the macroscopic quantum tunneling rate is required in order to simulate 
and understand the experiments. We present a derivation, based on the instanton technique, of an 
analytical expression of the escape rate in general quartic and symmetric sextic potentials compris- 
ing two escape paths. These new potentials were recently realized when creating a noise- insensitive 
phase qubit in the camel-back potential of a dc SQUID. 
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I. INTRODUCTION 

Due to the enormous progress in nanotechnology, su- 
perconducting nanocircuits are nowadays commonly ma- 
nipulated in the quantum regime^^^. The collective be- 
havior of superconducting electrons on the macroscopic 
scale is described with the superconducting phase of the 
Cooper pair condensate. As a function of this collective 
degree of freedom, the relevant potential energy for the 
phase is composed of minima separated by potential bar- 
riers. At sufficiently low temperatures, the quantization 
of energy in the potential wells generates quantum states 
that can be used in quantum information processing. The 
quantum dynamics as well as the readout of the quantum 
state are then governed by the tunneling between differ- 
ent minima of the potential. It is thus crucial to know the 
appropriate macroscopic quantum tunneling (MQT) rate 
of the phase for an arbitrary potential shape. Although 
this quantity can be calculated numerically in a general 
one- or multi-dimensional potential, it is much more con- 
venient to have an analytical solution when possible. In 
the following, we show that for a general quartic poten- 
tial and a symmetric sextic potential it is still possible 
to derive an analytical solution for the MQT rate. We 
will take a phase qubit as a physical example^, but our 
results apply to any system presenting such potentials. 

Usually, the precise shape of the potential depends on 
certain bias parameters. Most of the experiments on su- 
perconducting circuits are performed with current-biased 
Josephson junctions. The relevant potential for the phase 
is then the so-called tilted washboard potential^. If the 
bias current is close to the critical current, the wash- 
board potential is well approximated by a cubic poly- 
nomial with a local potential minimum from which the 
phase can escape into a continuum by tunneling. How- 
ever, interesting features appear away from this working 
point. In case of a SQUID loop for instance, at small cur- 
rent bias and with a flux bias close to half a flux quantum, 
the quantum states become weakly sensitive to current 



noise. In this region, the relevant potential has a quartic 
contribution and presents two escape points. This work- 
ing point was used in the experiment of Ref. 5, where a 
phase qubit was created in the so-called camel-back po- 
tential of a dc SQUID. For this experiment, the results of 
the present paper can be used for simulations and data 
processing^. 

In this paper we derive the expression of the MQT 
rate in a general quartic potential from which tunneling 
occurs through two possible escape paths, as depicted 
in Fig. [2j We use the instanton formalism developed in 
Ref. [7|. We also give the MQT rate for a symmetric sex- 
tic potential. The analytical expression is Eq. (|3T]) and 
the total escape rate turns out to be the sum of the es- 
cape rates for tunneling through each barrier. We find a 
simple result in the case of the experimentally relevant 
camel-back potential, Eq. (|35]) , depending on the bottom 
well frequency and the barrier height. We finally discuss 
the detailed shape of the potential in terms of the experi- 
mental parameters of the phase qubit in the two limiting 
cases of the washboard and the camel-back potential in 
Appendix [A] 



II. QUARTIC POTENTIAL IN A SQUID 

The SQUID is a superconducting loop interrupted by 
two Josephson junctions^, see Fig. [TJa. The Josephson 
junctions are characterized by their critical currents 7^ 2 
or equivalently their Josephson energies Ej^ .^ = ^/ci,2 5 
their capacitances Ci,2, and the phase differences (pi^2 
across them. Each arm of the loop has a self-inductance 
Li^2 and carries a partial current /i^2- The SQUID is 
biased by the current bias = Ii -\- 12 and the flux bias 
^5. We introduce the reduced variables x = ^{(fi + ^2) 
and y = ^{(fi — ^2)- In the following we consider a 
SQUID for which the capacitances in each arm are al- 
most equal, Ci €2- The Hamiltonian in the (x^y) 
plane reads H = ^{Ql + Ql) + Eju{x,y), where Qx 
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FIG. 1. (a) Schematic representation of the SQUID, (b) Po- 
tential of the SQUID u{x, y) and the path of minimum cur- 
vature connecting the minima through the saddle points, (c) 
Additional potentials \hjO± and \Mijo\Xq in units of the bar- 
rier height along the path of minimum curvature, the constant 
contribution has been removed, (d) Overlap of the ground 
state wave functions in the transverse direction as a function 
of the two abscissas zi and Z2. 



and Qy are the charges conjugate to the phases x and 
7/, = + is the total inverse capacitance, 

Ej = Ej-^ + Ej^ is the total Josephson energy. The char- 
acteristic frequency of the phase dynamics is the plasma 
frequency ujp = y/8EjEc /h^ where Ec = e^/2C is the 
charging energy of the SQUID. Due to the symmetry in 
the capacitances there is no coupling between the charges 
Qx and Qy. We define the reduced biases s = h/Ic^ 
yj) = $5/2^0 with the reduced flux quantum ^0 = ^/2e, 
the critical current asymmetry a = (/c2 ~ ^cj/^c, the 
loop inductance asymmetry r] = {L2 — and the 

junction-to-loop inductance ratio b = 2^o/^^c, where 
= Ici +^C2 is the total critical current and L = Li + L2 
is the total inductance. The reduced potential then reads 

u{x^y) = — cos X cosy— a sin X sin y—sx—sr]y-\-b{y— ys)^ ' 

(1) 

The dynamics is then equivalent to a fictitious particle 
of mass M = ^qC and coordinates (x, y) moving in the 
potential u{x^ y). 

The typical shape of the potential u{x^y), Eq. ([1]), is 
composed of local minima and saddle points connected by 
narrow valleys (see Fig. [TJb). The presence of extremal 
points is governed by the bias fields. Indeed, if s is larger 
than unity for example, there are no minima. For a given 
value of the fiux, we define the critical current as the value 
Sc at which the minima merge with the saddle points, 
consequently there are neither minima nor saddle points 
for bias currents s larger than Sc. The diagram sdys) is 
called the critical diagram. 

To study the phase dynamics from minimum to mini- 
mum in this potential, we will suppose that the fictitious 



particle follows the path of minimal curvature. We can 
then describe the dynamics with one variable, since the 
path in the plane (x, y) is determined by x and y = ({x) 
(see Fig. [TJb). We call z{x) the curvilinear abscissa along 
the path of minimum curvature, defined by 



z{x) = zo + / dxVl + C'(^0- (2) 

It is now tempting to consider the motion according to 
the effective one-dimensional potential Eju{z). In the 
following we check that the perpendicular dynamics does 
not generate additional terms. 

Let us consider a point (x, y) at abscissa z on this path 
and call the angle between the x axis and the direction 
of minimal curvature {6 G [— 7r,7r]) 

0{x) = - arctan — — — — . (3) 

2 yd^^u{x,y)-d^yu{x,y) J 

The direction of minimal curvature is denoted by X// and 
the perpendicular direction, given by the angle ^ + 7r/2, 
is denoted by X±. We can then express the Hamiltonian 
in terms of X//, X^, P// = MX//, and Pi_ = MX±_ by 
performing a rotation of angle 0. 

In the transverse direction, the frequency uj± = 

^Jd'^^u{x, y) Up is much larger than the bottom well 

frequency cuq in the direction of minimum curvature 
{oj±/oJo 7 in the experiment of Ref. 5). The per- 
pendicular dynamics is thus much faster than the dy- 
namics along the path of minimum curvature, and can 
be averaged out. The dynamics in the longitudinal di- 
rection is then governed by the effective Hamiltonian 
H// = {H)± obtained after averaging H over the trans- 
verse degree of freedom, the kinetic part being obtained 
using -\- Py = Pji + P^. To proceed, we consider an 
adiabatic evolution in the longitudinal direction and use 
X// as a fixed parameter. We suppose the perpendicu- 
lar dynamics frozen in the ground state of an effective 
potential of bottom well frequency uo^. The effective 
perpendicular dynamics, obtained within the harmonic 
approximation, is governed by the Hamiltonian 

H^ = ^Pl + \Mul{X^-Xof -\Mu:lXl (4) 

where the parameter Xq originates from the slope of 
the potential in the transverse direction. Averaging 
over the motion in the transverse direction gives rise 
to two contributions to the potential of JY//: {H±)± = 
^huj± — ^Muj\Xq. In practice, however, this additional 
potential {H±)±/Ej is only a small correction to u{z) 
(see Fig. [TJc) and will be neglected in the following. The 
quantum dynamics of the phase can thus be treated in 
the one- dimensional potential U{z) = Eju{z) following 
the path of minimum curvature. 

To quantify the adiabaticity, we evaluate the over- 
lap of the transverse ground-state wave functions at two 
points zi and Z2 along the path of minimal curvature. It 
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i)/2U4. 



(7) 



The point at an is always defined whereas ctl has a phys- 
ical meaning only when [/(Min^) < t/(Minc). When 
the coefficient U4, is larger than us, the potential is com- 
posed of two barriers. Because of these two "humps" (see 
Fig. [2]), this potential has been called the camel-back po- 
tential^. 

The dynamics in a given minimum is governed by the 
Hamiltonian 



aho 



(8) 



FIG. 2. General quartic potential U{z) — ^MuqZ^ (1 — U3Z — 
U4Z^), comprising two escape paths at an and ctl, passing 
through the barriers with a height Vl,r at the abscissa qL,R. 
The inset is the inverted potential obtained after a Wick ro- 
tation. 

turns out to be equal to R(zi,Z2) = v^N^C^i)^^/^^)]^' . 
If we choose Z2 = z -\- 5z close to zi = we have 
R{z^z + 5z) 1 — • order to obtain an 

adiabatic evolution in the ground state, the condition is 

thus ^^^^^ ^ 5z^ where the characteristic distance 5z is 

of the order of unity {Sz < 27r). The numerical evaluation 
of these quantities shows that the adiabaticity condition 
is satisfied (see Fig. [Hd). 

We choose the origin of the curvilinear abscissa along 
the path of minimum curvature at a particular minimum 
Mine with coordinates (xm, Vm)- Along the path of min- 
imum curvature {xX{x)), this minimum is connected to 
two saddle points SaddleL,i?, which are both connected 
to another minimum Miul^r. The shape of the barriers 
is given by the one-dimensional potential U{z). The con- 
vention is to call Miui^ the peripheral minimum with the 
lowest potential energy and to orient z towards Min^^. If 
we focus on the dynamics around the minimum, this po- 
tential is well approximated by its Taylor expansion up 
to fourth order, i.e. U{z) = U{0) + U{z) + 0{z^) with 

U{z) = \MujIz^ (1 - u^z - U4Z'^) , (5) 

where the bottom well frequency Uq ~ Up and the co- 
efficients ii3 > 0, 1^4 > are obtained numerically or 
with approximate expressions in specific cases (see Ap- 
pendix [Aj. The abscissa ql^r = 2;(SaddleL,i?) of the sad- 
dle points are 



qL,R = - 3ii3 T Y 9^3 + 32^/4 j /Su^. (6) 
We also define the escape points at abscissa aL,R where 



where Z = Muq /hz is the reduced position opera- 
tor, Pz is th e corresponding momentum operator, a = 
1^3 Y^a/Ma;o/2, and b = U4h/2MuJo. Due to the presence 
of the terms proportional to a and 6, Hamiltonian (j8j) 
describes an anharmonic oscillator. At sufficiently low 
temperatures below the plasma frequency, the energy 
spectrum of the particle is quantized in the minimum. 
Assuming the anharmonicity to be weak, a second or- 
der perturbation theory in a and b leads to a transition 
frequency between the levels n and n — 1 

hVn,n-l = ^^0 (1 - nKn) , (9) 

where A^ is the anharmonicity of the oscillator 

An = 36 + f a' + \b'{l7n + 7/n). (10) 

The anharmonicity depends on the working point 
{Ih^^h)' A sufficiently large anharmonicity is necessary 
to reach the two- level limit. Indeed, the manipulation of 
the states is performed with a microwave in resonance 
with the energy difference between the ground state and 
the first excited state and for a harmonic oscillator, where 
the anharmonicity vanishes, this microwave would excite 
all the levels. For a sufficiently large A, the two states |0) 
and |1) are decoupled from the others and constitute a 
qubit. The confining potential depending on the phases, 
this two level system is called a phase qubit 

To understand the quantum dynamics of the phase 
qubit in the camel-back potential, comprising two bar- 
riers, it is necessary to calculate the tunneling rate in a 
general quartic potential Eq. (j5j). We present the deriva- 
tion of the escape rate in the following section using the 
instanton technique. We also extend the result in the case 
of a sextic contribution when the potential is symmetric. 



III. DERIVATION OF THE MACROSCOPIC 
QUANTUM TUNNELING RATE 

The expression of the escape rate in the quartic po- 
tential sketched in Fig. [2l with a double escape path, is 
not explicitly given in the literature of MQT (see e.g. 
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Refs. Il3l - ll5[ ) and has to be calculated. We use the in- 
stanton technique^i^ to derive the escape rate from the 
metastable state of the quartic potential Eq. (|5]). 

The escape rate T is obtained from the quantum 
ground state of the fictitious particle confined in the well. 
Indeed, if Eq is the ground state energy, the norm of the 
particle wave function ^(^, t) in the well follows the tem- 
poral evolution 



21niEot/h 



which decays exponentially with the decay rate 



ImEo- 



(11) 



(12) 



The existence of a non-zero imaginary part comes from 
the possibility to tunnel out of the trap. The life- 
time of the state is equal to T~^. The quantum dy- 
namics of the particle is governed by the Hamiltonian 
H = ^P'^ + where the potential U{z) comprises 
a local but not global minimum, fixed at 2: = for con- 
venience. If we were dealing with a classical particle, the 
equilibrium state would be the particle at rest at the local 
minimum. This is however not possible in the presence of 
quantum fiuct nations. Moreover, quantum mechanically, 
the extension of the wave function allows the particle to 
evolve out of the trap where the potential energy is lower. 
The classical equilibrium is thus a "false ground state" 
The energy of the particle in the well is calculated from 
the propagator 



Git) = 



-iHt/h 



|0)=e 



-iEot/h 



(13) 



where |0) is the position eigenstate at z = 0. 

The probability amplitude Eq. (fT3)) is obtained after 
integrating over all possible trajectories z{t) satisfying 
the boundary conditions z{ti) = z{tf) = 0. For a given 
path, the phase of the contribution is the corresponding 
action in units of the action quantum H^^ . To proceed, 
let us perform first a Wick rotation: time t is replace by 
the imaginary time r = it. Then the action for a fixed 
trajectory becomes the Euclidean action 



S{z) 



T/2 



dr 



■T/2 



(14) 



where the dot means time derivative with respect to r 
and itij = /2. The effect of the Wick rotation on the 
action is thus to invert the potential (see Fig. [21 inset). 
In other words, in the part of space between the mini- 
mum and the exit point, the momentum of the particle 
is imaginary: = —2MU{z) < 0. But if we choose the 
time variable t as an imaginary variable r = it, then the 
motion is possible in the sense of classical dynamics in 
the inverted potential —U{z). In terms of path integrals 
the propagator reads 

rzi+T/2)=0 

G{T)=N / Vze-^^'^^^. (15) 



The parameter A/" is a normalization factor and Vz de- 
notes the integration over all functions z{t) obeying the 
boundary conditions. The propagator (p!5]) can be evalu- 
ated in the semiclassical limit, where the functional inte- 
gral is dominated by the stationary point z of S{z). The 
stationary point satisfies 5S{z)/5z\ - = 0, i.e., 



M't^U'{z) =0. 



(16) 



The action 

nT/2 

-M 



T/2 



dr 



of the stationary point is = 
' + U{z)\ . The path integral can be eval- 



uated with the method of steepest descent around the 
stationary trajectory, where the fiuctuations z can be ex- 
pressed on the eigenfunctions Zn of the operator — M9| + 
S,n^n' The functional integral, with 
Vz = Ylni'^^^)~^^in^ gives rise to a functional deter- 
minant. 



2(+T/2)=0 



Vze~^^^^^^ = e-^oT/2-So/h 



(-T/2)=0 



{det[- Md^^U'^z)]) (17) 



The prefactor JV is determined using the Gelfand-Yaglom 
formula, establishing that the ratio of two determinants 
of the same particle with the same energy in two differ- 
ent potentials is equal to the ratio of the corresponding 
wave functions. Using the asymptotic expression of the 
wave function of a harmonic oscillator of mass M and 
frequency ujq in the ground state, ipoiT/2) ~ e^^^ /2ujo, 
we get 



ttH 



-oJoT/2-So/h 



' det[- M ^U'\z)' 
det[-Md^^Muj^] 



(18) 
(19) 



Let us specify the shape of the trajectory z{r), solution 
of Eq. (p!6]) . to determine the action and the function 
K(T). We consider the nontrivial solutions where the 
particle can start at the top of the hill in the inverted 
potential, bounces off the potential on the right at 2; = 
(jR or on the left z = (Jl and returns to the top of 
the hill. For T ^ 00, we call this trajectory the right, 
and respectively the left, "bounce"^. The bounce has an 
energy £^0=0, thus 



and the bounce action reads 



5o,i = 2 / dz^2MU{z), 
Jo 



(20) 



(21) 



where the index j stands for L,R. To evaluate the func- 
tional determinant in Kj (T) , we follow the general frame- 
work of Ref. 7. The lowest eigenvalues have to be treated 



5 



22L] 

0.75- 
0.5- 

Tcb' 

0.25- 







U4 = 1 — Us 



FIG. 3. (Color online) Tunneling rate F, in red, for different 
barrier shapes, plotted in violet, from cubic to purely quartic. 
The rate is in units of the bottom well frequency cjo/Stt and 
the mass parameter is fixed to Mujo/h = 7. The tunneling 
rates Eq. (|33p for 1^4 = and Eq. (|35p for U3 = correspond 
to the dashed blue lines. 



separately. First, noticing that M z = U"{z)z^ we find 
that Z2 = ^s/M/Sojz is an eigenvector of the operator 
—Md^ + U"{z) with a vanishing eigenvalue A2. This 
eigenvalue is excluded from the determinant, now noted 
with a prime, and the integration over ^2 yields a prefac- 
tor ^JSQj/2iThMT. Second, as Z2 has a node, A2 is not 
the lowest eigenvalue. There is a negative eigen valu e, Ai, 
and the Gaussian integral over ^1 leads to i/2^\\i\. The 
intermediate expression for Kj{T) is then 



27rhM 



det' [-Md^^U'^z)] 



det[-MS2 + Ma;2)] 



T. (22) 



The negative eigenvalue guarantees a finite lifetime or 
equivalently a non-zero escape rate and characterizes 
the possibility for the particle to tunnel. The function 
Kj{T) can be evaluated from the asymptotic shape of 
the bounce orbit. If we suppose z(T ^ 00) ~ 
where 7^ = e^^ is obtained from 



Jo 



dz 



MujI 



1 



2U{z) 

the Gelfand-Yaglom formula gives finally rise to 



(23) 



(24) 



During the time T, several bounces occur, and the to- 
tal escape rate is obtained after summing over all the 
possible configurations. The multibounce configurations 
consist of riR right bounces and n l left bounces with cen- 
ters spread between — T/2 and T/2. We use the dilute 



instanton gas approximation, where the bounces are con- 
sidered to be independent, which is valid when the time 
between two successive attempts is larger than the tun- 
neling time. The action of a multibounce is ulSq ^tirSq 
and the contribution from the quantum fluctuations is 
K^^ /riLlriRl. The sum over the bounce number re- 
sults in 



G(T) = 



Mujo 
+e 



exp 



-UJqT 



The total decay rate F is 



Kr{T) . (25) 



(26) 



where 



F^- = 2e-'^°'^/^Imi^^-(T)/T. 



(27) 



As a conclusion, in the limit of a dilute gas of instantons, 
the total escape rate is simply the sum of the tunneling 
rates for each barrier. This simple result is a direct con- 
sequence of the dilute instanton gas approximation since 
interference effects between successive tunneling events 
would change the tunneling rate in one barrier and the 
resulting rate for two escape paths. 

As a result, the tunneling rate of a particle of mass 
M from an unstable state through the barrier potential 
U{z) with a double escape path reads F = 
where 





= Aje 




ll 












^0 



dz^j2MU{z). 



(28) 
(29) 
(30) 



Applied to the quartic potential of Eq. (|5]), this general 
result gives rise to the total tunneling rates 



TcrL,R 




X exp 



• (31) 



Eqs. (|26|) and (|3T]) constitute the central result of our 
work. The MQT rate is plotted in Fig. [3] for a potential 
varying from a cubic to a purely quartic profile. Start- 
ing from a purely cubic potential, the tunneling rate de- 
creases as the barrier height increases. When the tunnel- 
ing rates of both barriers are non- vanishing, the total tun- 
neling rate is enhanced and reaches another maximum for 
the purely quartic potential. This value of the tunneling 
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FIG. 4. (Color online) Effect of a sextic contribution on the 
MQT rate of a quartic potential. The rate is in units of the 
bottom well frequency uo/27t and the mass parameter is fixed 
to Mujo/h 7. 




FIG. 5. (Color online) Bounce orbits, Eqs. (|34| and {|36|, and 
MQT rates, Eqs. (|33)) and (|35)) . of the washboard potential 
and the camel-back potential. The tunneling time is of the 
order of the period To in the bottom well and, for the same 
bottom well frequency and barrier height, the escape rate is 
larger for the camel-back potential than for the washboard 
potential. 



rate is also obtained for a quartic potential when the cu- 
bic and quartic contributions are similar {us 0.5). 

When the potential is symmetric and has a sextic con- 
tribution, U{z) = \MujqZ^{1 — u/^z^ — uqz^), the ficti- 
tious particle can use two paths to escape at the exit 
points icTg with the same MQT rate rss/2. The exit 

points are defined with Gs = \J {y^u^ + 4:Uq — u/^)/2uq 

and the sextic coefficient satisfies uq > —u\/A {u^ > 0). 
At Uq = —2x4/4, the three minima are at the same height 
and a macroscopic quantum coherence develops between 
them. The possibility to tunnel between them lifts the 
degeneracy and the spectrum becomes composed of three 
hybridized levels separated by an energy given by the 
tunneling strength. Our general results can be applied 
to this symmetric sextic potential as well and leads to 



IV. WASHBOARD AND CAMEL-BACK 
POTENTIALS 

Applied to the washboard potential (2x4 = 0), the stan- 
dard MQT tunnehng rate is recovered^^ 



wb 



^ 27r V hjQ 



exp 



36 V 

5 hujQ 



(33) 



where V = 2Mcj^/27c|. From Eq. (|20]), the bounce orbit 
is found to be 



^wb(T) 



0'L,R 



cosh (cjot/2) 



(34) 



In the case of the camel-back potential {us = 0), we 
get the tunneling rate 



=4cJr 



Mujo 



7Th ^2^ 

2^X4 



ul - 



4ii6 /2 
■ arctan 



U4 



}■ 

(32) 



In Appendix [A] we show that the camel-back potential 
is in fact better fitted with a sextic contribution. The 
resulting MQT rate is presented in Fig. |4l The sextic 
contribution mainly sharpens the camel-back potential 
and hence increases the tunneling rate. 

The tunneling rate of the ground state can be used to 
calculate the tunneling rate for the excited statesi^. 



cb 



16 



2^ 



/27rV 


16 


V ' 









(35) 



where V = Mco^/Sci. Eq. ^ can also be obtained from 
the symmetric sextic MQT rate Eq. (|32]) with uq = 0. 
The bounce orbit of the camel-back potential is equal to 



0'L,R 



cosh(cJoT) 



(36) 



As seen in Sec. [Ill these two regimes can be observed 
in a phase qubit close to the critical line and around 
the working point biased with a weak current and half 
a flux quantum, respectively. The bounce orbits as well 
as the MQT rates are plotted in Fig. O The tunneling 
time turns out to be of the order of the period of oscil- 
lation To = 27rujQ^. The approximation of a dilute gas 
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of instantons is valid when the time between two tun- 
neling events is larger than the tunneling time Tq, 
i.e. r < cJo/27r. This corresponds to a barrier height 
larger than the ground state energy Jiwo/2. The escape 
rate from the camel-back potential of parameters (cjo, V) 
corresponds to the escape rate from a cubic potential of 
parameters (^cjq ^ 0.89a;o, OMV). For two 

potentials of the same frequency and same barrier height, 
obtained when C3 = y/cj, the ratio of the escape rates 

A 28 V 1 n V 

is Fcb/Fwb = 575^1^^-0 - 0.77e ^^-0. The tunnel- 
ing rate out of the camel-back potential is then larger 
than the escape rate from the washboard potential when 
V/huJo > ^ log ^ 0.14, see Fig. El 

Finally, for a purely sextic potential U{z) = 
^MuJqZ^{1 — uqz^)^ the total rate reads 



I2ttV 


Sv^TT V ' 




4 hbjQ 



(37) 



where the barrier height is equal to V = MuOq/'^^/^uq for 
Uq > 0. 



V. CONCLUSION 




x/tt 



FIG. 6. (Color online) Potential (top panel) along the path of 
minimum curvature y = ys — cos(x)/26 (bottom panel), for 
a dc SQUID with parameters a = 0.0072, ry = 0.69, b = 3.03, 
and the biases s — —0.0062 and ys — 0.53l7r. Top panel: 
quart ic and symmetric sextic expansions, in red and green 
respectively. Bottom panel: potential u{x,y) and position 
obtained from numerics of the minima and saddle points, red 
and green crosses, respectively. 



Using the instanton technique we calculate the escape 
rate from the metastable state of a quartic and a sym- 
metric sextic potential, in which the fictitious particle 
can tunnel through two different barrier potentials. This 
work is of interest for Josephson junction based nanocir- 
cuits, such as in the field of circuit quantum electro- 
dynamics, where the potential shapes are designed for 
an efficient quantum information processing. In particu- 
lar, the results can be directly used in the case of phase 
qubits, where we give the tunneling rate as a function 
of the biases in two different regimes, namely close to 
the critical line and also close to the working point with 
a weak current bias and half a ffux quantum magnetic 
bias. 
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Appendix A: Topography of the phase qubit 

In this Appendix we present the shape of the one- 
dimensional potential, i.e. cjq, 1^3, and 1^4, of a phase 
qubit in two different regimes, namely the washboard po- 
tential and the camel-back potential. The bottom well 
frequency and the barrier height can be directly used to 



obtain the MQT rate as a function of the experimental 
parameters. 



1. Washboard potential 

Close to the critical line the one-dimensional potential 
has a washboard shape. We Taylor expand the potential 
in X, y and s at third order around Xc, yc^ and Sc in 
the direction of vanishing curvature, given by the angle 
tan^c = —dxx^{^c^yc)/d'^yU{xc^yc)- The potential can 
be expressed in terms of the barrier height V and the 
bottom well frequency ujq as follows 

Uiz) =\mu:Iz- (1 - y^z) , (Al) 

(A2) 

/ \ 1/4 

WO = (^ j (cos^e + r/sin^e)'/'(l - sf'^UJp, 

(A3) 

where C3 = cos (9c (1 + 2sin2|9c)(sin Xc COS yc — 

acosxc sin?/c) + sin^c(l + 2 cos^ 6>c) (cos Xc sin — 
asinxc cos?/c). These expressions generalize the known 
results for a symmetric SQUID^^ (j] = 0) and can be 
directly used in Eq. ([33|) to obtain the MQT rate. 
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2. Camel-back potential 



where 



We focus on the domain of bias fields around the 
critical point {sc = —a^ysc = 7r/2 + 1/26, = 
0,?/c = 7r/2}, the others can be obtained by symme- 
try. The path of minimum curvature follows the tra- 
jectory y = — cos{x)/2b. The Tailor expansion of the 
one-dimensional potential up to fourth order in z reads 
U{z) = \MujIz'^(1 - usz - U4z'^) with 



Us 



3C2 ' 



U4 



12C2' 



(A4) 



C2 



1 + 7^5 — COS P 



2b 



sin /3, 



cs = a cos P 



3a 



C4 



2b 

4 + 775 — 7 cos P 
2b 



sin /3, 



(A5) 
(A6) 
(A7) 



and P = ys — VBc- The path of minimum curvature and 
the corresponding potential with a quartic and a sym- 
metric sextic expansions are plotted in Fig. [6l The sextic 
contribution is necessary to fit accurately the potential 
between the two exit points. 
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